Possible pairing symmetries in SrPtAs with a local lack of inversion center 
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We discuss possible pairing symmetries in the hexagonal pnictide superconductor SrPtAs. The 
local lack of inversion symmetry of the two distinct conducting layers in the unit cell results in 
a special spin-orbit coupling with a staggered structure. We classify the pairing symmetry by the 
global crystal point group D^d, and suggest some candidates for the stable state using a tight-binding 
model with an in-plane, density-density type pairing interaction. We may have some unconventional 
states like s + /-wave and a mixture of chiral d-wave and chiral p-wave. The spin orbit coupling 
is larger than the interlayer hopping, and the mixing between spin-singlet and triplet states can be 
seen in spite of the fact that the system has a global inversion center. 

PACS numbers: 74.20.Rp 
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The relation between crystal structure and pairing 
symmetry plays an important role in unconventional 
superconductivity^ Pairing states can be categorized 
with respect to the irreducible representations of the 
point group of the crystal lattice and do not mix unless 
they belong to the same representation. Since the Pauli 
principle requires that the momentum part of singlet and 
triplet states possess even and odd parity, respectively, 
their mixing is prohibited in a system with inversion 
symmetry. Superconductivity in non-centrosymmetric 
systems, i.e., CePtsSi, opens however the possibility of 
singlet-triplet mixing^— It plays a key role to explain the 
puzzling behavior of the observed nuclear spin-lattice re- 
laxation rate £ Microscopically, this mixing is caused 
by an anti-symmetric spin-orbit coupling (SOC). 

Recently, possible singlet-triplet mixing in centrosym- 
metric systems with a local lack of inversion symmetry, 
such as special crystal lattices or heterostructures, was 
discussed^ The recently-discovered hexagonal pnictide 
superconductor SrPtAs 8 (T c = 2AK) belongs to the for- 
mer case of a special crystal structure. The unit cell pos- 
sess a global inversion center and its point group is D^d- 
There are two distinct honeycomb Pt-As layers within 
the unit cell each of which has no inversion center. LDA 
calculations revealed that these two layers are conduct- 
ing with only a small inter-layer hopping, i.e., the system 
is quasi- two-dimensional (quasi-2D). In addition, a large 
splitting of the bands due to anti-symmetric spin-orbit 
coupling (SOC) was seen. The consequences of this lo- 
cal lack of inversion symmetry on magnetic properties 
of the superconducting phased as well as on electronic 
phenomenal has previously been studied. In this work, 
we aim at clarifying its role for the pairing symmetry. 

Table Q] shows the classification of the pairing states 
based on the global symmetry of the crystal D^- Wc 
assume intra-layer pairing due to the quasi-2D nature of 
the system, and focus on on-site and nearest-neighbor- 
site (nn-site) pairing interactions. It is intriguing that 
in this table both even-parity spin-triplet and odd-parity 
spin-singlet pairing appear. The reason is that we have 



two distinct layers in the unit cell indicated by I = 1,2, 
and we can introduce an odd-parity factor (—1)' under 
the global inversion operation. Multiplying this factor to 
a certain pair wave function results in even-parity spin- 
triplet or odd-parity spin-singlet states. Moreover, spin- 
singlet and triplet states coexist in some irreducible rep- 
resentations, namely A\ g , E g , Am and E u . Therefore, 
mixing of spin-singlet and triplet states becomes possible 
in these representations despite the parity conservation. 

Since there is no experimental information on the pair- 
ing symmetry at present, we discuss some potential can- 
didates for the stable symmetry within a simple model. 
We use a tight-binding description for electrons on the 
Pt sites with a Hamiltonian consisting of two parts: 
H = Hq + H sc . The first part, H , is the one-body 
Hamiltonian introduced by Refs. [H and [if] in order to 
reproduce the LDA band structure of SrPtAs, 
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where c kls (c kl J ) is the annihilation (creation) oper- 
ator of an electron in the 6-th band (b = 1,2,3) 
with crystal momentum k, spin s in the l-th layer 
(/ = 1,2). In the above equation, we introduced 
<7° (f°) and a 1 (f 1 ), the unit and Pauli matrices act- 
ing on the spin (layer) space. Including Pt nearest- 
neighbor hopping within the plane, as well as nearcst- 
and next-nearest-neighbor hopping between the planes, 



one finds e 
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4 J2n cos ^ ' ^« + ^c2 cos(ck z ), and 
4 fc) cos(^c/2)[l + exp(-ifc • T 3 ) + exp(ifc • T 2 )\ 
with Ti = (0,o,0), T 2 = {V3a/2, -a/2,0), and T 3 = 
(— V3o/2, — a/2, 0) the in-plane nearest-neighbor bond 
vectors used in the tight-binding approach (a and c are 
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TABLE I: (a) Spin-singlet, and (b) spin-triplet basis gap functions. This classification is based on D^d symmetry. The 
index I = 1,2 denotes two distinct layers. The definitions for functions of crystal momentum k are, e*, = ~^2 n cosk ■ T n , 
e fc = ]C n ^™cosfc • T n , o k = J] n sinfe • T n , o+ = ^ n a;™sinfc • T n , e~ - e+*,o~ = o+*, where T n= i >2 ,3 is the bond vector 
between nearest-neighbor sites, and cj n = exp[2n7ri/3]. Note that we have even-parity spin-triplet part and odd-parity spin- 
singlet part due to the odd-parity factor ( — 
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in-plane and inter- layer lattice constants). An important 
ingredient is the locally anti-symmetric SOC Xki, which 
reads A& = z J2 n sin k ■ T n for each band. This term 
is symmetric under global inversion, but anti-symmetric 
under the local inversion operation in each layer. Due 
to the Kramers degeneracy, there are only two branches 
in the energy spectrum of the Hamiltonian (JlJ for each 
band 
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We use tight-binding parameters from Ref. which lead 
to Fermi surfaces as shown in Fig. Q] With this param- 
eters, the outermost band, labelled band 3, is the domi- 
nant band with 74% of the total density of states (DOS) 
due to its proximity to the van Hove singularity (vHS) at 
the M points in the Brillouin zone (BZ). Note that the 
ratio Oib/tc , which parametrizes the effect of the local 
lack of inversion symmetry, is comparable or larger than 
1. This large ratio plays an essential role for the mixing 
between spin-singlet and spin-triplet state, as we will sec 
below. 



For the pairing term H sc in the total Hamiltonian we 
assume intra-layer interactions including density-density 
type attractive interaction, as well as inter-band pair 
scatterings allowed by the kinematics. Using the basis 
functions from Table U H sc is written in Fourier form as 

H -Vl/ bM ' b ' <k k')r im r (h)t c (W) c (b ' ] (A\ 




FIG. 1: Fermi surfaces at (a) k z — 0, and (b) k z — n/c. 
Inner blue, middle red, and outer green lines show the Fermi 
surfaces for band 1, band 2, and band 3, respectively. Note 
that there are two branches in each band as suggested in Eq. 
©, but one of the branches in band 3 does not cross the 
Fermi level at k z =0. 



with 

r 

bb' ,(r,m)* /« . » \ fr\ 

ikl jk'l' \ O ~i O ~y)siS2\°~y ~j)S3S4-i \P) 



where and g^ n are the coupling constants for on-site 
and nearest-neighbor channels. The pairing instability in 
this model occurs in band 3 with its dominant contribu- 
tion to the DOS. Smaller gaps then open on the other 
two bands due to pair scattering. 

We solve the linearized gap equation (the eigenvalue 
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equation for T c 



where the sum runs over repeated indices, and 
G^(k,iuj n ) = {a g> (io;„f + e^ b) ) + a fc A fc • 6- ® f 3 } V) 

is the normal-state Matsubara Green's function. All the 
possible gap functions arc listed as 
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where A^ and A r ^ m=1 2 are the order parameters, and 

s^p and tp are the mixing ratios of subdominant spin- 
singlet and triplet parts, respectively. We see in T = A\ g 
and A2u that there is a mixing between on-site and 
nearest-neighbor-site pairings, besides the spin-singlet 
and triplet mixing. We neglect the band dependence of 
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and introduce repulsive inter-band in- 
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teractions g^ {nn) = g Q ' n{nn) = -0.05, keeping <? on(nn) 
This choice is motivated by the nesting-like structures 
between band 2 and 3, and band 1 and 3, respectively^ 
We can then calculate the state with the maximum eigen- 
value T™ ax at a point (g n, g nn ) in the coupling constant 
space. 

Figure [5] shows the obtained phase diagram. The A\ g 
state is stabilized when the on-site attraction is dom- 
inant, whereas the Am state becomes stable in the pa- 
rameter region where the nn-site attraction is comparable 
to, or larger than the on-site coupling. From Table U and 



FIG. 2: The phase diagram of the stable pairing states in the 
coupling constant space (<?on,gnn). The tight-binding param- 
eters suggested by LDA calculation^^ is used. The sequences 
of dots show equal T c lines at T c /t[ 2) = 10~ 5 , 10~ 4 , 10~ 3 from 
bottom to top. 



both Hebel-Slichter peak and power-law type tempera- 
ture dependence of T{~ like CcPt 3 Si£ The gap structure 
involves sign changes which give rise to zero-energy An- 
dreev bound states at certain surfaces, e.g. for the normal 
vector [010] ~ Note that the relation of the bound state 
and t opolog y of the wave function has been discussed in 
Refs. Il3ljl4l . This state belongs to the class AIII of the 
topological classification-^. 

The locally anti-symmetric SOC introduces a mixing 
between spin-singlet and triplet parts, which is propor- 
tional to 
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have "s + /" -wave pairing symmetry with the s-wave (/- 
wave) component dominant while the /-wave (s-wave) 
component with an odd-parity factor (—1)' is subdom- 
inant. Therefore, the quasiparticle excitations are fully 
gapped in the A\ g state, whereas line nodes appear in 
the Am state. The A 2u state invokes a full coherence 
factor due to the s-wave component, and would show 



This suggests that the mixing is suppressed by a large 
inter-layer hopping, as expected, since the system has 
global inversion symmetry and the locally anti-symmetric 
nature is smeared out when the three dimensionality 
becomes strong. Such a behavior is also seen in the 
magnetic susceptibility^ In this system, however, the 
inter-layer hopping has been estimated to be compara- 
ble or smaller than the SOCS^ and we hence expect 
a finite value of mixing. Indeed, around the bound- 
ary between the A\ g and Am phases in Fig. [2] we find 
enhanced mixing ratios. Their magnitudes are almost 
band-independent and typical values are (sj£ , ) = 

z(b) l9 Jb) 19 . _ 



we see that both, the A lq and the A 2u state, (-0.51,0.12) in the A lg phase, and (s^,s 



(0.15,-0.18) in the Am phase (definitions of the ratios 
are given in Eq. ©). 

Figure [3] shows the phase diagram for a shifted chemi- 
cal potential such that band 3 approaches the vHS. The 
enhanced DOS naturally leads to reduced coupling con- 
stants for the same T c as compared to the previous sit- 
uation. More remarkably, the E g state shows up in the 
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FIG. 3: The phase diagram of the stable paring state in the 
coupling constant space (<?on,(/nn) in the DOS enhanced sit- 
uation, where the Fermi level is located at the vHS point 
of band 3. The sequences of dots show equal T c lines at 
T c /if ) = 10" 5 , 10" 4 , 10" 3 from bottom to top. 

region where the on-site coupling is repulsive. One of the 
reasons for its stability is that the amplitude of the sin- 
glet component l^fY™! nas peaks at the saddle points, 
which is compatible with the Fermi surface structure. 
This phase involves two degenerate basis states indicated 
by to = 1, 2 in Eq. ([5|). and they make up a Kramers pair. 
A fourth-order analysis of the Ginzburg-Landau theory 
yields to degenelate states (A^, ^A^' 2 ) = (1, 0), (0, 1), 
which both break time-reversal symmetry. We focus here 
on the first configuration and set A^/ 2 = 0- Expand- 
ing the spin-singlet component around the zone-central 
axes k x = k y = gives V'f'z' 1 = (k x + ik y ) 2 with 
d x 2_ y 2 +id xy -w&ve symmetry, or chiral d-wave symmetry. 
Note that d^_„2 and d X y components arc degenerated in 
the three-fold rotational symmetry. The same expansion 
for the spin-triplet part yields d^j" 1 = (— l) l (k x — ik y )z 
with chiral p-wave symmetry like S^RuO^i 1 ^ The chi- 
ral d-wave part has L z = +2, whereas the chiral p- 
wave part L z = — 1 (L z : z-component of the relative 
angular momentum of the pair). These states can mix 



with each other as indicated by Table |Tpi The mixed 
state is classified into class A in the scheme of the topo- 
logical classification^ Due to the chiral nature of the 
pairing, this state has a non-zero value for the Chern 
number defined by the vorticity of the quasiparticle wave 
function in k spac o 17 ' 18 and supports chiral edge states 
topologically i 19 ! 20 

Our analysis provides insight into the basic trends of 
the hexagonal superconductor SrPtAs whose electrons 
experience a locally non-centrosymmetric environment. 
The A\ g state is stable in the electron-phonon coupling 
limit where on-site attraction is dominant . On the other 
hand, in the strongly-correlated limit with on-site re- 
pulsion or strong nearest-neighbor attraction, the A2 U 
state is stabilized. In this state, line nodes coming from 
the spin-triplet component cause a power-law behavior of 
Tj~ , whereas a Hebel-Slichter peak arises slightly below 
T c due to the coherence factor of the spin-singlet compo- 
nent, in analogy with CePt3Si£ Such a behavior would 
be a strong signal of the locally anti-symmetric SOC. As 
mentioned, the nodal structure results in Andreev bound 
states at a certain surface^ which is related to the topol- 
ogy of the bulk state d 3 ' 14 The E g state is possible in 
some particular cases like DOS enhanced situation ow- 
ing to the vHS of the saddle points in the hexagonal BZ. 
This state breaks time-reversal symmetry whose signal 
could be detected by /xSR measurement for spontaneous 
magnetization around impurities and also the Kerr rota- 
tion experiment, for examples. The state has chirality 
which is characterized by the Chern number, and leads 
to topologically-protected chiral edge statesj 19 i 20 
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